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Abstract 

Cohesive energies for twelve cubic III-V semiconductors with zincblende 
structure have been determined using an ab-initio scheme. Correlation contri- 
butions, in particular, have been evaluated using the coupled-cluster approach 
with single and double excitations (CCSD). This was done by means of incre- 
ments obtained for localized bond orbitals and for pairs and triples of such 
bonds. Combining these results with corresponding Hartree-Fock data, we 

recover about 92 % of the experimental cohesive energies. 
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I. INTRODUCTION 



Ab-initio calculations for ground state properties of solids are most often performed 
nowadays within the density-functional (DFT) scheme. Typically the exchange and correla- 
tion contributions are calculated within the local-density approximation (LDA)§, or, more 
recently, within the generalized-gradient approximation (GGA)i. These methods yield quite 
good results for cohesive energies and other properties, but do not provide many-body wave 
functions; thus, a systematic improvement towards the exact result appears to be difficult. 
Hartree-Fock self-consistent-field (HF-SCF) calculations for solids, on the other hand, while 
lacking electron correlation as per its definition, have the merit of treating non-local exchange 
exactly and of supplying a good starting point for subsequent correlation calculations. One 
possibility suggested in literature is to calculate correlation contributions on top of HF 
within LDA or GGaI, but the potential advantage through the availability of a many-body 
wavefunction is absent and the results do not seem to achieve a significant improvement over 
exchange-correlation (xc) -DFT . 

For finite systems like atoms and molecules, accurate quantum-chemical methods such as 
configuration-interaction (CI) or coupled-cluster (CC) approaches have been developed for 
the determination of correlation energies. Infinite systems such as solids require the use of 
size consistent approximations for treating electron correlations like coupled electron pair 
approximations (CEPA) or CC approaches. Because of the local character of the correlation 
hole one can expand the correlation energy in terms of local incrementsi. The idea is thereby 
to determine the required matrix elements by studying local exitations in clusters which are 
accessible to a full quantum-chemical treatment. A formal derivation of the method for an 
infinite periodic system can be found in Ref. [5] or [6] when the CEPA scheme is used. 
The derivation for the CC method is very similar, see, e.g., Ref. [7]. The method of in- 
crements has previously been applied to ground-state properties (cohesive energies, lattice 
constants and bulk moduli) of elementary semiconductors and has led to very good and 
reliable resultsi. 
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In the present paper we extend the method to polar semiconductors, especially to the twelve 
cubic III-V compounds. In order to assess the validity of the calculated correlation energies, 
we evaluate their differential contributions to the cohesive energies of these materials. For 
a comparison with experiment, however, we also need reliable SCF ground state energies. 
In Section II we perform the SCF calculations using the program package CRYSTAL92i. In 
Section III we sketch the method of correlation-energy increments and report on computa- 
tional details. The results are presented and discussed in Section IV. Conclusions follow in 
Section V. 



II. HARTREE-FOCK CALCULATIONS 

For discussing the influence of electron correlation effects in solids, reliable SCF calcu- 
lations are prerequisites. For the compounds which we are interested in, SCF calculations 
with the program package Crystal88 have been performed by the Torino groupi. However 
small split-valence basis sets have been used in these calculations, and it is not clear how 
closely the results approach the HF limit. 

Here, we want to supplement their results with our own calculations employing the scalar- 
relativistic energy-consistent pseudopotentials (PP) of the Stuttgart group together with 
the corresponding atomic basis sets0. The cores replaced by the pseudopotentials are 
and for group 13 and 15, respectively, with two exceptions: 3-valence-electron pseu- 
dopotentials for the post-d elements Gallium and Indium have been found to underestimate 
the closed-shell repulsion of the underlying d shell on valence electrons of neighbouring 
atomsBEl; we therefore performed the SCF calculations for Ga and In compounds with 
13-valence electron pseudopotentialsS, explicitly treating the highest occupied d shell. 
For the solids, we generated {4s4pld)/[3s3pld] Gaussian valence basis sets (Ga: from 
the atomic basis set0 {8s6p7d) the very diffuse exponents are neglected in obtaining 
{7s4:p6d)/[3s3p2d], In: from the atomic basis set0 {6s6p6d) the very diffuse exponents are 
neglected in obtaining (4s4p5(i)/[3s3p2(i]) as follows: starting from the energy-optimized 
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atomic basis sets derived for the PP just mentioned00, the inner two (Ga: five) functions 
of s and p symmetry were contracted using atomic ground-state orbital coefficients. The 
outer two s and p functions were left uncontracted. The most diffuse exponents of the atom 
have little effect in the solid, where, due to close-packing, their role is taken over by basis 
functions of the neighboring atoms. Moreover, the use of basis functions that are too diffuse 
leads to numerical problems in Crystal. Therefore we re-optimized the two outer s and p 
functions as well as an additional d polarization function for the solid. The crystal-optimized 
basis sets are listed in Table |l[ 

Using the above basis sets, we determined SCF ground state energies of the solid with 
CRYSTAL92i. When evaluating cohesive energies, we subtracted corresponding energies 
of free ground-state atoms obtained with the original atom-optimized set sS§ using the 
quantum-chemical ab-initio program MOLPRO940. The results for the SCF cohesive ener- 
gies at experimentally determined values of lattice constants^ are listed in Table and 
comparison is made to experimental cohesive energies corrected by the phonon zero-point 
energies |kB6D (derived from the Debye modefl) as well as by atomic spin-orbit splittingsEl. 
It is seen that binding energies, at the HF level, are between 50% and 70% of the exper- 
imental values leaving room for significant correlation contributions. The agreement with 
the previously published SCF cohesive energieJi is very good for compounds containing B 
and Al, while small deviations of ~2% of the cohesive energy are found for Ga and In com- 
pounds. The stability of the results with quite different pseudopotentials and basis sets 
(Causa et al used PP's of the Toulouse group and [2s2pl(i] basis sets§) is gratifying to note, 
and we expect to be not very far from the basis-set limit of the SCF cohesive energy with 
our [3s3pl(i] valence basis sets. 

III. CORRELATION EFFECTS 

Calculating correlation effects using an expansion in terms of local increments^ is formally 
similar to treating the hierarchy of nth order atomic Bethe-Goldstone equations^. Here we 
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only want to sketch the basic ideas and some important formulae (for more details see Ref. 
[5] and [6]). The method relies on localized bond orbitals generated in a SCF reference 
calculation. One-bond correlation-energy increments are obtained by correlating each of 
the localized orbitals separately while keeping the other ones inactive. In the present work 
we are using the coupled-cluster approach with single and double substitutions (CCSD). 
This yields a first approximation of the correlation energy 

-^corr = Xl^«' (1) 

i 

which corresponds to the correlation energy of independent bonds. 

In the next step we include the correlations of pairs of bonds. Only the non-additive part 
Aejj of the two-bond correlation energy e^j is needed. 

Acy = - (ci + ej). (2) 

Higher order increments are defined analogously. For the three-bond increment, for example, 
one has 

Aeijfc = €-ijk - (ci + + Cfe) - (Aejj + Ae^fc + Aeife). (3) 

The correlation energy of the solid is finally obtained by adding up all the increments with 
appropriate weight factors: 

E^}^ = E + ^ E Ae^. + \ E A^^.^ + •••• (4) 

i ij ijk 

It is obvious that by calculating higher and higher increments the exact correlation energy 
within CCSD is determined. 

The procedure described above is only useful if the incremental expansion is well convergent, 
i.e. if increments up to, say, three-bond increments are sufficient, and if increments decrease 
rapidly with increasing distance between localized orbitals. These conditions were shown to 
be well met in the case of elementary semiconductorsi, but have to be checked again for the 
polar compounds. 



Since (dynamical) correlation is a local effect, the increments should be fairly local entities 
at least for semiconductors and insulators. We use this property to calculate the correlation- 
energy increments in finite clusters. We select the clusters as fragments of the zincblende 
structure so that we can calculate all two-bond increments up to third nearest neighbors 
and all nearest-neighbor three-bond increments. The clusters are shown in Fig.|I[ The bond 
length between the group 13 (X) and the group 15 atom (Y) is taken to be the same as in 
the solid at the experimental equilibrium lattice constant. The dangling bonds are saturated 
with hydrogens. The X — H and the Y — H distances, respectively, were optimized in CCSD 
calculations for XYHe clusters, yielding dBK= 1.187 A, dj<su= 1.016 A, dAm= 1.614 A, dpu= 
1.408 A, dGa.H= 1.621 A, dAsR= 1.525 A; for In and Sb, we adopted the SnH4 distance (1.711 
A). Instead of constructing localized bond functions from the Bloch states derived using the 
Crystal92 program, it is simpler to determine them from a cluster calculation. One merely 
has to show that, for a given increment, convergency is reached when the size of the cluster 
is sufficiently large. With this in mind we performed standard SCF calculations for each 
cluster and localized the bonds according to the Foster-Boys criterion^ within the occupied 
valence space in Ci symmetry. Following the procedure described above we calculated the 
correlation-energy increments at the CCSD level (using the program package Molpro940) 
successively correlating more and more of the localized X — Y bonds. 

For hydrogen we chose Dunning's0 double-zeta basis without a polarization function. Two 
different basis sets are used for the other elements. Basis A is of the same quality as for 
the SCF calculations of the preceding section: (4s4p)/[3 s3p]E valence basis sets energy- 
optimized for the atoms, with d polarization functions optimized in CCSD calculations for 
XYHg clusters (see Table H). Extended basis sets (B) have been generated by uncontracting 
the sp functions of basis A and by replacing the single ci-function by a 2dlf polarization 
set. The latter was obtained by adding a second more diffuse d function to the one of basis 
A and by optimizing the / exponent in CCSD calculations for the XYHg clusters (see Table 

0)- 

We have checked the convergence of the incremental expansion for all polar substances and 
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will discuss it here for GaAs. As can be seen from Table |ITT| , the increments decrease quite 
rapidly with increasing distance of the bonds. The nearest neighbor two-bond increments 
contribute 51% of the correlation energy of GaAs, all next-nearest neighbors 13% and the 
third-nearest ones only 6% . After these increments we truncate the expansion because the 
fourth-nearest neighbors' contribution is only 0.5% . Note that the integrated truncation 
error decreases with r~^, because the interaction is van der Waals like and the number 
of pairs of bonds grows like r^. 

The convergence with respect to the number of bonds correlated simultaneously in the 
incremental expansion is also quite satisfactory. For GaAs again, the one-bond increment 
contributes 38% of the correlation energy, the two-bond increments are very important 
(70%) but yield a correlation energy too large in magnitude, while the three-bond increments 
reduce it by 8%. An estimate of the nearest-neighbor four-bond increments and next-nearest 
neighbor three-bond increments is ~0.5% . Thus, the overall error due to truncations in the 
incremental expansion is less than 1% of the correlation energy, the same accuracy as was 
reached for the elementary semiconductorJl. 

The next point we have to check is the transferability of the increments from clusters to the 
solid. We test the sensitivity of the chemical surroundings by determining the increments 
in different clusters. The change from GaAsHg, where the Ga — As bond is surrounded only 
by hydrogens, to Ga4As4Hi8, where the inner Ga — As bond has the same neighbors as in 
the solid, is 0.005 a.u. for the one-bond increment. The differences in the nearest-neighbor 
two-bond increments are -0.005 a.u., so that the main errors due to lack of transferability 
cancel. The tendency for individual (one- and two-bond) errors to be of different sign is 
seen in all III-V compounds. We estimate an upper limit for the transferability error to 1% 
of the correlation energy. 

Thus, we end up with a total error of 2% for approximations inherent in our use of the 
incremental expansion of the correlation energy. Not included are shortcomings of one- and 
many-particle basis sets in the determination of individual increments; the largest part of 
this remaining error is probably due to limitations of the one-particle basis set and will be 



discussed in the next section. 



IV. RESULTS AND DISCUSSION 

Applying the method of increments as described in the preceding section, we determined 
correlation contributions to cohesive energies for all twelve cubic III-V semiconductors. The 
increments were always taken from the largest possible cluster (cf. Table III and Fig. ^ 
and multiplied by the weight factors appropriate for the zincblende structure. The corre- 
lation contributions to the cohesive energies were obtained as -Eco" = -^soM ~ -^atom i 
per unit cell. The results for the two different basis sets are shown in Table For basis 
set A we obtain an average of 86% of the experimental cohesive energies, which amounts 
to ~67% of the 'experimental' correlation contributions to the cohesive energies (defined 
here as the differences between the experimental cohesive energies and the corresponding 
SCF values). The larger basis set B yields a substantial improvement, ~82 % of the cor- 
relation contributions to the cohesive energies. Overall we obtain with basis B an average 
of 92% of the cohesive energy with only small fluctuations for the different materials. In 
earlier calculations for the elementary semiconductorsi, we obtained ~94% of the cohesive 
energy. The deviation may be due, in part, to the different correlation treatment. For the 
elementary semiconductors we used a coupled-electron-pair approximation (CEPA-0). It is 
well known that CEPA-0 yields larger correlation contributions as the CCSD approach; in 
test calculations we found a deviation of ~1%. If the basis set, especially the polarisation 
set is enlarged even more, e.g., to 3d2flg, a further increase of the correlation effects can 
be expected. In test calculations for GaAs, where the one-bond increment and the nearest 
neighbor two-bond increments were determined with this extended basis set, we obtained 
an increase of 10% of the correlation contribution to the cohesive energy. 
For comparison, we have also listed in Table ^ results from the literature which have been 
obtained with other methods. LDaS overestimates the cohesive energies by ~20 % , GGA 
is at the same level of accuracy or slightly better than our results. In contrast to these 
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methods, our approach lends itself to systematic improvement; moreover, a detailed insight 
into the nature of the correlation effects becomes possible. 

V. CONCLUSIONS 

We have determined cohesive energies of twelve cubic III-V semiconductors, both at the 
SCF and the CCSD levels. The SCF results have been obtained with the Crystal92 code 
using relativistic energy-consistent pseudopotentials and [3s3pl(i] Gaussian valence basis 
sets. Electron correlations are described at the CCSD level with the method of local incre- 
ments, using the same pseudopotentials in cluster calculations with the ab-initio program 
MoLPRo94. Increasing cluster and basis-set size, this approach allows for a systematic 
improvement of accuracy towards the fully correlated solid-state limit. The results show 
that the method works well for polar III-V semiconductors with zincblende structure: the 
calculated cohesive energies are 92±3% of the experimental values using an uncontracted 
(4s4p2dl/) valence basis set. Work is underway in our laboratory to apply the method to 
other ground-state properties of the III-V compounds. 
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TABLES 

TABLE L Crystal-optimized basis sets for the IILV semiconductors. Of the s and p functions, 
only the two outer uncontracted exponents are listed. The inner contracted ones are chosen to be 
the same as in the original atom-optimized basis setJi^'0. 
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The very diffuse d-projector of the boron pseudopotential was neglected. 



TABLE II. Polarization functions used in the basis sets for the CCSD calculations. 
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TABLE III. Correlation-energy increments for GaAs (in a.u.), determined at the CCSD level 
using basis set A. For the numbering of the clusters and bonds involved, see Fig. || . 
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TABLE IV. Cohesive energies per unit cell (in a.u.), at different theoretical levels (cf. text). De- 
viations from experimental values (in percent) are given in parentheses. For comparison, literature 
data from DFT calculations are also reported. 
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FIGURES 

FIG. 1. X„Y„HTO-clusters treated at the CCSD level. Big numbers designate clusters, small 
numbers the bonds in each cluster. H-atoms are not drawn. 
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